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1. INTRODUCTION 


The Agriculture and Resources Inventory Surveys Through Aerospace Remote 
Sensing (AgRISTARS) is a 6-year program of research, development, evaluation, 
and application of aerospace remote sensing for monitoring agricultural 
resources beginning in fiscal year (FY) 1980, The AgRISTARS program is a co- 
operative effort of the National Aeronautics and Space Administration (NASA), 
the U.S. Agency for International Development (AID), and the U.S. Departments 
of Agriculture, Commerce, arH the Interior (USDA, USDC, and USDI). 

The goal of the program is to determine the usefulness, cost, and extent to 
which aerospace remote sensing data can be integrated into existing or future 
USDA systems to improve the objectivity, reliability, timeliness, and adequacy 
of information required to carry out USDA missions. The overall approach com- 
prises a balanced program of remote sensing research, development, and testing 
which involves domestic resource management as well as commodity production 
information needs. 

The technical program is structured into eight major projects as follows: 

1. Early Warning/Crop Condition Assessment (EW/CCA) 

2. Foreign Commodity Production Forecasting (FCPF) 

3. Yield Model Development (YMD) 

4. Supporting Research (SR) 

5. Soil Moisture (SM) 

6. Domestic Crops and Land Cover (DC/LC) 

7. Renewable Resources Inventory (RRI) 

8. Conservation and Pollution (C/P) 

The majority of these projects will make direct use of information on crop 
phenology. Phenological information will be pertinent to these several areas 
of the projects including classification, acreage and yield estimation, and 
detection of episodal events. 
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Recent research at a number of locations has been directed toward developing 
techniques for correcting variations in Landsat spectral signals caused by 
local differences in growth stages. The general technique consists of fitting 
ari exponential equation to multitemporal spectral data. One result of such 
fits is an estimate of the date on which the maximum spectral response would 
have been observed if there had been daily spectral observations. 

The original algorithms of this type were developed by Badhwar (1979, ref. 1). 
Work conducted at the Environmental Research Institute of Michigan (ERIM) by 
Cicone et al . (1979, ref. 2) resulted in the use of a different functional 
form for the fitted equation. For crop stage estimation, a simplified 
algorithm was developed. 

This paper presents the results of a preliminary study of the characteristics 
of the ERIM algorithm. A methodology for using the calculated date of peak 
response to estimate growth stage was used to evaluate the algorithm and to 
provide a comparison with the Badhwar approach. 

2. BACKGROUND 

The ERIM crop calendar shift algorithm operates on the greenness component of 
the Kauth-Thomas transformation (Kauth et al., 1976, ref. 3). The day of 
peak response is determined by examining the cross-correlation between the 
observed values and the corresponding values of a reference profile. The pro- 
gram requires that three spectral acquisitions separated by at least 15 days 
in the emergence to preharvest period be available. The existing algorithm 
checks for these conditions and does not require manual selection of acquisi- 
tions. The data should be preprocessed to remove the effects of satellite 
calibration, Sun angle, and atmospheric haze. ERIM recommends the use of its 
spatially varying XSTAR algorithm and associated correction factors for 
preprocessing. 
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The basic algorithm was provided in the form of a Fortran subroutine (Horvath, 
1980, ref. 4). Complete documentation of the program along with a user's 
guide and a series of test cases for program verification is given by Crist 
and Malila (1980, rnf. 8). 

For this study, a data set was constructed using spectral data collected over 
spring wheat test sites during the 1978 crop year in the U.S. northern Great 
Plains. Periodic ground-truth observations were available for selected fields 
on actual growth stages. Table 1 shows a summary of the data set. From 1 to 
13 fields were available in each of 22 locations. A total of 681 growth stage 
observations were available. This data set was used previously (Cate et al., 
1980, ref. 6) to evaluate the Badhwar crop calendar shift algorithm. Manual 
quality assurance had already been performed on the data. As a result, the 
only preprocessing employed in this evaluation was satellite calibration and 
Sun-angle correction. 

3. ESTIMATION OF PEAK SPECTRAL RESPONSE DATE 

The ERIM algorithm provides estimates of the date of peak spectral response. 
The program provides the value of the cross-correlation coefficient for those 
fields which have satisfactory acquisition histories. Table 2 summarizes the 
results of initial model runs. Estimates were made for approximately 65 per- 
cent of the fields in the original data set. The number of fields, number of 
Landsat acquisitions, and the average and range of observed R 2 values for the 
cross-correlation associated with the estimated peak response date are shown. 
Most fields had R 2 values in excess of 0.9, indicating that the profile fits 
were of a high quality. No geographical pattern was found in the distribution 
of the goodness of fit as measured by the R 2 values. 

4. ESTIMATION OF DEVELOPMENT STAGE 

Since the basic algorithm only provides estimates of the date of peak spectral* 
response, an attempt was made to develop a technique for the estimation of 
individual growth stages. The ground- truth observations were taken on the 
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TARLC 1.- SUMMARY Of DATA 
THE TRIM CROP CALENDAR 
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TABU 2.- SUMMARY OF MODEL PERFORMANCES BY FIELD 


Segment 

Number of 
fields 

Number of 
acquisitions 

; Cros^corrrlatlon* R 2 


Range 

1392 

5 

3 

0.944 

0.873-0.995 

1394 

4 

4-5 

.978 

.974- .984 

1457 

5 

3-4 

.982 

.938- .999 

1461 

11 

5-6 

.935 

.916- .963 

1518 

11 

4 

.858 

.599- .990 

1537 

1 

5 

,970 

.970 

1542 

6 

3-5 

.932 

.824- .984 

1544 

2 

4 

.863 

.752- .973 

1553 

2 

5 

.798 

.753- .842 

1566 

2 

3 

.748 

.649- .846 

1584 

4 

3 

.986 

.980- .990 

1636 

y 

3-5 

.940 

.870- .980 

1668 

2 

3 

.992 

.992 

1811 

10 

3-4 

.911 

.757- .980 

1825 

13 

4 

.979 

.932- .998 

1920 

5 

3-4 

.919 

.894- .944 

1924 

2 

5 

.967 

.948- .985 


I 


I 

j 

j 


i 

i 
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Feekes scale (Lame, 1^41, ref, 3 1. In previews wo*h , c^'p et al. {inar% 
ref. H) developef' a rpiatSontM; he 4 wpp’" ,* he ‘-^ek-ps *.fnlr Arid the pprront, of 
the normal growing * rsnn w H oh Had nlanspd at each s:agn. T ahle 3 gives o 
recent revision of thi *■ relationship {Kodqrs, iQRP, ref, f.l. This tab’e 
provides a vehicle for 1 ransforndno the ‘ppkrs ohspr*'*t ivs info a quantify 
which is linear in time, The table show* the percent n* the plant growth 
Which had elapsed at the midpoint of path foekes staue well as the staop 
duration. 

The estimate of day of peak oreenness was calculated 4 *or each available field. 
A regression betwoer the ground-truth stage in percent ana the days from date 
of peak greenness was run for the entire data sr:. r fc-‘ir r T ar ana’ vs is was 
carried out using peak date estimates obtained from me Badhwar aloer^ihm. 
Table 4 gives the regression results, eiots for each algorithm are shown in 
figures 1 and 2. 

The estimated growth stage for a qiven field on a particular acquisition was 
calculated using the following equation. 

Stage (?) o a r + a 1 { n KDAY - AC.mftY) 


where 

Stage (%) * percent of the growing season 

ap.a’L * intercept, and slope, respectively, the regression 

PKDAY « date of peak spectral response for the field 
ACQDAY -< acquisition date of the field ground truth 

Table 4 indicates that the models a«e estimating substantially different dates 
of peak spectral response. The FRF algorithm peaks when HB percent of the 
growing season has elapsed as opposed to fi?, percent for the Badhwar algorithm. 
Table 3 shows that the peaks of the two algorithms are roughly a week apart; 
the peak of the ERF algorithm occurs near the beginning of heading, and the 
peak of the Badhwar algorithm occurs near the end of heading. 
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TABLE 3.- RELATIONSHIP BETWEEN FEEKES DEVELOPMENT STAGE 
AND PERCFNT OF NORMAL SFASON GROWTH FOP SPRING WHEAT 
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Figure 1.- Plot of development stage versus days from time of peak spectral response 

for the ERM crop calendar shift algorithm. 





Figure 2.- Plot of development sta*. versus days from time of peak spectral response 

for the Badhwar crop calendar shift algorithm. 





During preliminary review of these results, it has been suggested that this 
variation may have been the result of differences in the data preprocessing 
employed with the two algorithms. Spectral greenness values often are nega- 
tive during part of the growing season. The usual procedure is to add a 
constant value to greenness tn ensure positive numbers. A constant value of 
32 has been adopted. The ERIM algorithm subtracts 25 from the greenness, 
resulting in an effective offset of 7. The Badhwar algorithm was run with an 
offset equal to 32. This difference may partly explain the discrepancy 
between the two models. 

5. EVALUATION OF DEVELOPMENT STAGE ESTIMATES 

The ERIM algorithm calculated peak response dates for 92 of the 139 fields 
listed in table 1. A total of 489 of the possible 681 growth stage estimates 
were made. The Badhwar algorithm estimates were summarized for the same 
field-date combinations. The ground truth was converted to percent of season 
and the error of each estimate was calculated in percent. Table 5 compares 
the results obtained from the two algorithms summarizing the bias and root 
mean square error (RMSE) by growth stage. In general, the algorithm's per- 
formances were similar with the Badhwar algorithm showing higher RMSE for most 
stages. 

The overall results were summarized in table 6 to provide a better measure of 
model performance. Two methods of summarization are shown. In the first 
case, the errors were weighted by the number of observations in each stage. 

All growth stages were used. The second method examined only the period 
between emergence and ripe. In this case, the errors were weighted by the 
duration of the stage as given in table 3. The ERIM alqorithm again shows the 
best results, particularly in the second case. Since spring wheat has approx- 
imately a 100-day growing season, these errors may be interpreted broadly as 
errors in days. Thus, the ERIM algorithm results in a 1- to 1.5-day reduction 
in error for estimating the date of a given growth stage when compared with 
the Badhwar algorithm. 
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TABLE 6.- SUMMARY OF RESULTS FOR THE ERIM AND BADHWAR 
CROP CALENDAR SHIFT ALGORITHMS WEIGHTED BY NUMBER 
OF OBSERVATIONS AND BY DURATION OF STAGE 


Algorithm 

Feekes 

stages 

Weighting 

technique 

Bias 

RMSE 

ERIM 

All 

Number of observations 

-0.7 

9.92 

Badhwar 

All 

Number of observations 

.5 

10.79 

ERIM 

1.0 to 11.4 

Stage duration 

.9 

8.63 

Badhwar 

1.0 to 11.4 

Stage duration 

2.3 

10.42 
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Additional runs were made to determine the effect of relaxing the acquisition 
restrictions in the software. When a 9-day separation in acquisitions was 

allowed, a bias of 0.3 and an RMSE of 10.22 was observed, and the number of 

estimates increased to 590. When the definition of the growing season was 
increased to 110 days, as opposed to 90 days in the basic algorithms, the bias 

was -0.1, the RMSE was 10.38, and the number of estimates was 613. When a 

110-day season and 9-day acquisitions were allowed, the bias was 0.3 days, the 
RMSE was 10.59 and the number of estimates was 681. Comparing these results 
to the baseline results — which showed a bias of -0.7, an RMSE of 9.92, and 
489 estimates — indicates that, while some deterioration in performance occurs 
when the restrictions are relaxed, the increase in the number of estimates 
probably justifies the loss in accuracy. 

6. RECOMMENDATIONS 

The ERIM algorithm is recommended for testing with independent 1979 crop year 
data. The model should be tested in both the original configuration and with 
the relaxed restrictions, allowing 9-day acquisitions and an expanded defini- 
tion of the growing season. 

The performance of models of this type can probably be enhanced by using 
meteorological data from the date of peak spectral response to estimate the 
dates of other stages. It is recommended that the use of accumulated degree 
days be investigated for this purpose. One simple approach along these lines 
would be to develop an equation analogous to the one presented in the method- 
ology section, which relates stage to the difference between the estimated 
peak day and the acquisition date. 
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